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Abstract 

We obtain conditions for the existence of an attractor in the system of equations describing a tachyon 
warm inflationary model with bulk viscosity taken into account. When these conditions are met the evolution 
approaches slow-roll regime. We present the primordial power spectrum for the tachyon field by considering 
a dissipation coefficient depending on the scalar field and temperature. 
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I. INTRODUCTION 


The inflationary paradigm has proved to be successful in solving the problems related with the 


initial conditions of the standard cosmological model (horizon, flatness, monopoles) 


it provides a natural explanation to the large scale structures we see today [2] and the anisotropy 


B. 


Q. 


Besides 


observed in the cosmic microwave background 

In the standard inflationary scenario a single scalar field dominates the dynamics driving a phase 
of accelerated expansion in the very early universe. During the inflationary period the universe 
results supercooled and a subsequent phase of reheating is necessary. An alternative scenario is 
denominated warm inflation [4]. In this model the dissipative effects are important during inflation 
and radiation production occurs concurrently with the expansion. The dissipative effect is mediated 


by a friction term describing the decaying of the scalar field into a therma’ 


bath and the fluctuations 
At the end of a period 


a: 


of the scalar field arise from thermal rather than quantum contributions 
of warm inflation the temperature of the universe is high enough to smoothly go into the radiation 
dominated phase of the standard cosmological model. 


Several models of warm inflation has been studied [6y, particularly a model of warm inflation 


driven by a tachyon field was proposed in [7]. A tachyon field [8| is an interesting candidate to 
drive inflationary expansion 0 and it is associated with unstable D-Branes in string theory [lo| . 

Prom the theoretical point of view, any inflationary model must be a dynamically plausible 
solution to the cosmological equations, i.e. a viable inflationary scenario must be an attractor in 
the phase space of solutions to the cosmological field equations fll]. 

Usually the slow-roll approximation is used in inflationary models because the equations are 
considerably simplified. Nevertheless, it is not always verified that the slow-roll solution is an 
attractor solution to the dynamical set of equations describing the system. 


B * 


In Refs. [12- 16] the consistency of some warm inflationary models has been revised. Particularly, 

m RefQ the cons,stency of a tachyon warm ™ a 
dissipation coefficient independent of temperature. In Ref. [171] the primordial spectrum of density 
fluctuations for a warm tachyon model was calculated taking into account a term of viscous pressure. 

Viscous pressure may significantly influence the dynamics of warm inflation and it is well mo¬ 
tivated on physical grounds. A term of viscous pressure may arise from the decay of heavy fields 
into light fields or by the production of particles by the inflaton 

The aim of this work is to show that the slow-roll solution for tachyon warm inflation with 
bulk viscous pressure is behaved as an attractor under particular conditions and to present the 


ist: in 

Q. 
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primordial power spectrum of scalar perturbations and the corresponding spectral index for a 
tachyon warm inflationary model with a dissipation coefficient depending on the temperature. 

The paper is organized as follows. In section [IT] we present a model of warm tachyon inflation 
with viscous pressure. In section HITI we perform a linear stability analysis, we look for the existence 
of an attractor in the space of parameters describing the model. In section [TV] we study the 
associated primordial power spectrum by considering a dissipation coefficient with temperature 
dependence. In section lYll we present some examples with an exponential potential for the tachyon 
field. Finally in section IVlII we summarize our findings. 


II. BASIC EQUATIONS 


We consider a spatially flat, homogeneous and isotropic universe with expansion rate H. The 
matter content is given by a tachyon field cj> and a radiation bath of entropy energy densityjs and 
temperature T. The dynamics of tachyon warm inflation with viscous pressure is given by 


/ _r c 


3 H 2 = 


1 


V 


Ml 


+ sT 


= — — \J\ — u 2 u — 3 Hu — 


Vj 


1 — u 2 


Ts = Yu 2 — 3H(Ts + II) 


(1) 

( 2 ) 

( 3 ) 


where we have conveniently defined the dimensionless quantity u = <j>. The evolution is governed 
by the effective potential V(4>, T) and a dissipation coefficient Y((j>, T) (T > 0 by the Second Law of 
Thermodynamics). We consider the bulk viscous pressure given by the relation II = —3 (H, where 
the phenomenological coefficient of bulk viscosity £ is assumed as a function of the energy density 
p and it is positive defined. The subscript in a function denotes derivative, we work in natural 
units and we consider 8nG = Mp 2 . In Eqs. dT]{3]) , the tachyon field has dimension [</>] = M _1 , the 
dissipation coefficient [T] = M 5 , and the bulk viscous pressure [II] = M 4 , where M denotes the 
dimension of mass. 

The energy density of the tachyon field is defined as p $ = . We note that the weak energy 

condition for this fluid impose u 2 < 1. On the other hand, the energy density of the radiation 
bath is related to the entropy density by r yp = sT where 7 corresponds to the adiabatic coefficient, 
which must slowly vary with expansion but for the sake of simplicity, we will consider it constant 
and in the range 1 < 7 < 2. 

From the thermodynamic relation s = — ^ , where / is the Helmholtz free energy, we can relate 
the entropy density with the effective potential as s ~ -4^ Q. 
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Basically, the slow-roll approximation consists of consider a dominant potential and neglect 
terms of the highest order in time derivatives from Eqs.dlD-Q. The validity of this approximation 
will depend on a set of dimensionless slow-roll parameters defined in terms of the potential as: 


e = 


M 2 V 2 


_ 71 ^2 

r\ = M P - 


P = Mj 


V*Y, 


(4) 


2 V 3 ’ '' V 2 ’ r ±r p V 2 T 

In order to give account of the temperature dependence in the effective potential and in the 


dissipation coefficient it is necessary to introduce two additional dimensionless parameters 


Q: 


b = 


4>T 

V<h 


c = 


TtT 


(5) 


During the inflationary period the Hubble expansion rate is nearly constant | jp | <C 1, this 
defines the first slow-roll parameter ep which we can write from Eqs. dH)-({3]) as: 


eH = 


H 

H 2 


3 u 2 + ( 


sT i n 

V V 


)(wT 


- u 2 - 1 - 


Vu z 

3 HV 


1 - v / E r 


2(1 + fVT^ 


( 6 ) 


The second slow-roll parameter is given by: 

H V 


VT - 


u* 


VH = 


+ 


HH HV Hu 1 - u 


1 - 


3 HV 


sT 

V 


3 u 2 + + £) ^ 4\/1 - u 


— u 2 — 1 I — 


r u 2 

3 HV 


1 + VvH 



(2« 2 (-3 + - i))i + -1) - 1>) 


3« 2 + (f + S) (4vT 


- u 2 - 1 ) - 4^7 (1 - a/ 1^72 


The smallness of these parameters give the relative size of terms neglected in the slow-roll 
approximation. We note that in order to have ep 1 (or accelerated expansion) we must impose 
u 2 -C 1, sT <C V, n <C V and Yu 2 <C 3 HV] furthermore the condition 7 /# « 1 implies u -C Hu, 
T -C ET and n -C HE. Additionally we impose that the decay of tachyon held into radiation 
is quasi-stable, s <C sH. The end of inflation is presented when ep becomes 1 or accelerated 
expansion ends. 

In the slow-roll regime Eqs. ([T])- d3j) reduce to: 

V 


3 H z = 


M 2 




3H(1 + r)u + = 0 

3 H(sT + n) - Yu 2 = 0 


(7) 

( 8 ) 
(9) 


where r is a dimensionless coefficient defined as r = jyy- This coefficient allow us to distinguish 
between two regimes of warm inflation, the weak dissipation regime (r < 1 ) and the strong 
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dissipation regime (r 2 > 1). We will focus on this latter case because in the weak regime viscosity 
contributions are unimportant [18]. 

From Eq. ([9]) we note that a consistency condition for II < 0 is that | -y | < 1. 

It is worth to mention that in the slow-roll and strong dissipation regime Eqs.©-© look alike 
the slow-roll equations in the case of warm inflation with a canonical scalar field. 

The number of e-foldings is given by: 


1 V 2 

N (<t>) = -TTT / — (! + r )# 


Mf - . 

where the subscript e denotes the end of inflation. 


V d 


( 10 ) 


III. STABILITY ANALYSIS 

Following Ref. [t 3| . we perturb Eqs.| 
system in matrix form as, 


and (|3|) around the slow roll solution xq and express the 


8X' = M{xq) 8X - x' 0 (11) 

where for convenience we have considered derivatives with respect to </> denoted with primes and 

A B \ I Su 


M(x 0 ) = 


5X = 


and 


C D 


5s 


XQ 



( 12 ) 


The term x' Q is present because the slow-roll solution is not an exact solution to Eqs.([T])-([3|). 

In order that the slow-roll solution becomes an attractor to the system, the eigenvalues of M{x o) 
must be negative and the components of vector x' 0 must be small. The condition on the eigenvalues 
is met when: the determinant of matrix M(x o) is positive and the trace is negative. 

Given that the slow-roll solution is not an exact solution to Eqs. dH)-d3j) the term x' 0 in (fill) must 
be small. The size of x' 0 depends on the quantities, 


s 

3 t 

' / 2r + 8 

Hs 

A ' 

VV r + 1 

u 

3 1 

( f ( 3 “- 

Hu 

A 1 



-e - 27? + 

XJ r 


r ~ 1 a , 3 , (1 + 7") 2 3q (1 -r 2 ) 

-P + ^bxi— - ~i bc x^— 


(13) 


(1 + r) 


ii Xl i i 

+ ( co; - — ] r? + 


<prk + \ bc( i +r)x 'i (14) 


where we define cr = -4p. % = 1 + a, xi = 7(1 + ic), u = 7 — 1 and A = 3y(r + l)xi — 3w(c(l — r) + 


2 b(r + 1)) + yy?:^1 + yyy-xj • We note that the conditions |\ <C 1 and | \ <C 1 are satisfied 
in the strong dissipation regime provided that: 

|e| < r, |r?| <C r, \/3\ < r, | 6 | < 1 (15) 
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For tachyon warm inflation with viscous pressure we have: 


A = 
B = 


H 

u 

H 


—3(1 + r) + 


e(4r + 1) 


C = 5(6 + 6 a- 


D = 


u 

H 


(1 + r ) 2 

- 6(1 

e(l + 2 er) 
(1 + r ) 2 


n , . , . „ 3 _ re (7 — 1) f a(l — a) 4(7 — 1 ) 

— 3(7 — 1) (cr — 6(1 + ?’)) + -aaa + . , 2 I ---F 3c H 

s \ 2(1 + r) z \ 7 — 1 7 


u \3( 7 — l)( c ~ 1 + <rc) - 3(1 + 7 *1) - \p<r* + 5^5 (V/3 _ 2 ™< 7 + V 


( 16 ) 

(17) 

(18) 
(19) 


where we define a = 12 + 7 ^ 117 24 ^ , /3 = (2 37 ^ 2 7 ^ and a = S = v and we have considered 

C (p) oc P l . 

The requirement of positive determinant of matrix M in (11211 is more restrictive than the 
constraint over the trace, we find that in order to the slow-roll regime behaves as an attractor in 
the strong dissipation regime (r 1 ) we must impose the following condition: 

27(1 + at) — /3ad — 46(1 + <r )(7 — 1) 


c < 


( 20 ) 


2 (l + a)( 7 -l) 
where we have considered 6 > 0 , |cr| < 1 and |cr/| < 1 . 

We note that the parameter c is not required to be small and the condition on the parameter 6 
is more restrictive than the conditions over the other parameters. It is worth to mention that the 
condition on 6 restricts the temperature corrections of the potential to be small. 


We note that the result for a canonical scalar field 1121 is recovered in the limit cr —» 0, and the 
result in Ref. [131] for a canonical scalar field with viscous pressure is reproduced by taking l —>• 0. 


IV. DENSITY FLUCTUATIONS 

The authors of Ref. Q have determined the primordial power spectrum of density fluctuations 
for the case of a tachyon field with a term of bulk viscous pressure: 




H 2 (l +r ) 2 ^ 2 

where (Scj)) 2 are the corresponding fluctuations of the tachyon field and 

r (T/U), 

F W>) = - J [A{4>) + dcj), where A(<i>) = 3H + r /y > 

2H + T/V 




8 G(<p) (3H + r/vy 


r + ahv - 7 -1 + 


CA (In V) ( j ) \ (InU) 


C J 3 7 il ( 3 H + T/V) J V 


G(4>) = 1 - 


2 sT + 3 n + i (sT + n) (uf -1 


8 H 2 


( 21 ) 


( 22 ) 


(23) 


(24) 
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The fluctuations of the tachyon field are generated by thermal interaction with the radiation 


field, rather than quantum fluctuations. In the case of a dissipation coefficient depending on the 
temperature it is hard to get analytical solutions. The authors of Ref. 2o| give an approximate 
expression for their numerical result of the power spectrum for a canonical scalar field. These 
results are derived to leading order in the slow-roll approximation. It is possible to develop an 
analogous procedure in the case of tachyon warm inflation with viscous pressure. To leading order 
in slow-roll the perturbation equations for the tachyon field and the radiation bath can be simplified 
as it is shown in the Appendix. The contribution of the viscous pressure is parametrized with the 
a parameter defined in section um which we assume approximately constant during inflation. 


= 


-HTy/l + r ( 1 + 


rg 


(25) 


where <5 = ^+ 7 ia) an d fg is obtained numerically (see Appendix). 


rg 

1 = 0 

l = 0.5 

l = 1 

c = 1 

8.53 

13.81 

41.70 

c = 2 

7.66 

14.78 

33.93 

c = 3 

7.27 

15.85 

25.88 


TABLE I: Numerical values of rg, we assume 7 = 4/3 and a = —0.1 


The Eg. ( 1251) is consistent with the spectrum of the scalar field with a dissipation coefficient 

n □ 

independent of the temperature [5(. We note that for II = —3<Zo H the result in Ref. [20] is valid in 
the case of viscous pressure. The result in Eg. (1251) is valid for 5 > 0 or a > — 

From Eqs. (l2ip and (I25|) we calculate the new spectral index in the strong dissipation regime, 


n, — 1 = 


d In <5 


„ , = bpi + m + VP3 + PPi ~ 2 ( 8 - 2e + -e 


H 


where we define 8 =-,e =— e and 

37(1 + <r)(l + la) (4c o + 2clo + 3c 2 y + 2(2 + c)l^a) 


Pi = 


P 2 = 


P 3 = 


Pi = 


16(w + / 7 c) (7 + I'ya + cw(l + a)) 

3c 2 cuj(2 + 3a) — 2(cj + Z7<t)(57 — 2 + 3lja) — c (8w 2 — 97 s + a (6cu 2 — lj(8 + 7 ) + 


(26) 


(27) 




4(cu + I'ya )(j + l-ya + cw(l + <r)) 

3 c 2 <u 7(1 + a) — 2ccu(l + a)(ui + I'ya) — 4(c o + l^a){\ + (1 + (—1 + l)j)a) 


2 (ui + l^a)(yy + I'ya + cw( 1 + o - )) 
3c7 2 (1 + la) + 2 (a; + I'ya){2 - 7 + (2 + (-2 + l)^)a) 


(29) 

(30) 


4(cu + I'ya )(j + I'ya + cu(l + a)) 

We have considered that the interval of wave number is related with the number of e-foldings N 
by the relation din A: = dN. 


7 


















V. TENSOR PERTURBATIONS 


Tensor perturbations do not couple to the thermal background, so gravitational waves are only 

m , 


generated by quantum fluctuations, as in standard fluctuations 

2 V 


2 f H_ V _ 

^ n) “ 6 t r 2 M|, 


(31) 


The tensor spectral index is given as usual by ut = —2e. 

The tensor-scalar ratio for tachyon warm inflation in the strong dissipation regime is given by, 


where Vr = ™5 2 H . 


V T _ 2VZ IH2F 
Vr tt 2 ^ T 6 


(32) 


VI. EXAMPLE 

In order to maintain the validity of the slow-roll approximation as an attractor in the strong 
dissipation regime, the condition on the slow-roll parameter b guarantees that the thermal cor¬ 


rections to the effective potential must be small. This allow us to consider separable potential 
T ) = Vi(T) + V2(4>). On the other hand, by using Noether symmetry the authors of Ref. 


as 


22 ] 


found that if a tachyon field is driving accelerated expansion, then the potential associated to the 
field must be exponential. The above arguments motivate the use of the following potential: 


V{4>) = V 0 e 


= v nf ,~ a 4> 


(33) 


where Vo and a are positive constants and for the sake of simplicity we have only considered a 
potential depending on the scalar field. 

Taking into account the results of section HVl we consider a dissipation coefficient of the form: 


r (<f>,T) = r 0 e~l a ^T c 


(34) 


where To and c are constants positive defined. 

In this exmaple we also consider II = —£o pH. 

Usually, the end of a period of slow-roll inflation is implemented through the condition e# = 1 
or equivalently by imposing a null acceleration to the system. It is not always possible to find this 
behaviour because en must be an increasing function in order to reach the condition at the end of 
a long enough period of inflation [16]. In the case of a bulk viscosity coefficient proportional to the 









radiation density, it is possible to analytically show that a dissipation coefficient depending on the 
temperature allows to achieve this behaviour. From Eqs.Q-Q and Ea. (|34D we have: 


H 

CH = ~ip = 


Vo are 


2 PO ar t> 


c r r, 


r 1 0 


2 r n 


M 2 p V* ,2 a 2 


(y/SMpj _ 3e~^^/VoCo] 


4+c 


(35) 


where we assumed that the energy density is related to the temperature through the relation 
p = C r T 4 . Given that cf) is an increasing function, ep is an increasing function provided that 


- 3Co H 


a = — 


< 1 , i.e. €h < 1 during inflation and one at the end of inflation. 


From Eqs.([7|)-([9]) and considering p = C r T 4 we can explicitly find: 


T(H) = 


r(H) = 


f \/3C r To 
V MpVn /2 a 2 


1 

4+c 


-(7-3Co H) 


\M P V C 

MpTq f y/3C r T o ( \ 

VZVq 12 \M p Vq / 2 o . 2 7 0 ) 


4+c 


(36) 

(37) 


We note that during inflation the temperature is a decreasing function, so is p. 


Given that in our example we consider the temperature dependence in T, it is not possible to find 
more explicit results. Several examples where T depends only on cj) were studied in Ref.[l3|. The 
current constraints on the amplitude of the primordial power spectrum for scalar perturbations 
Vji(ko) = 2.23 x 10 ~ 9 and the tensor to scalar ratio R(ko) < 0.11 [23j] (for a pivot scale fco = 


0.002Mpc 4 ), allow us to constrain the value of the potential at the Hubble horizon crossing to be 

v; < io- 8 m 4 . 

In our example the slow-roll parameters are given by: 

.,2 


- - P 

e = r >=3 = 


and we can rewrite the spectral index in 

n s - 1 = 


M 2 a“ 

2 rV 


as: 


< 1 , 6 = 0 


Mna 2 


2 rV 


P 2 ~\ +P3 + 3+4 - 2) 


(38) 


(39) 


We note that n s — 1 is always negative for l = 1, 7 = 4/3 and c = 1,2,3, —0.25 < a < 0. 
Consequently in a scenario of tachyon warm inflation with viscous pressure it is common to have 
a red tilted primordial power spectrum for scalar perturbations. 

VII. FINAL REMARKS 


We analysed a model of tachyon warm inflation by taking into account a term of viscous pressure. 
In the strong dissipation regime we found conditions in order to the slow-roll solution behaves as 
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an attractor to the dynamical system of equations describing the model. The slow-roll regime 
is implemented by imposing restrictions on the size of the five slow-roll parameters of the model 


associated to the effective potential and the dissipation coefficient. The conditions on e, r/. (3 
and b are the standard ones for warm inflation [ 12 ], but the condition on the c parameter has 
contributions associated to the viscous pressure term and it is less restrictive than the case without 
viscous pressure [l(| for — 1 < er < 0 and l < 1 . 

We calculated an approximated expression for the primordial power spectrum of a tachyon 
field driving a period of warm inflation. We consider a temperature dependence in the dissipation 
coefficient and viscous pressure. This result follow from Ref.j^oj] and the numerical integration of 
the corresponding perturbed equations is expected to be performed in future work. 

We presented the primordial power spectrum for density fluctuations by considering a dissipation 
coefficient with temperature dependence, and we calculate the corresponding spectral index to this 
model. We work through a specific example with a decaying exponential potential for the tachyon 
field. We found that it is possible to build up a consistent model with temperature dependence in 
the dissipation coefficient, which properly finish the inflationary period. 
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Appendices 

For tachyon warm inflation with viscous pressure Eqs.(45) and (46) of Ref. [20] become: 

- (3r + 2)^ + 5cj) + = y/W^/T£ 

5p” — ( 37(1 + la — c/4) + 4)^- + (157(1 +la - c/4))^ + (u + 7 la)8p r + (T <fi)8(j) = 0 
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In calculating the primordial spectrum of the tachyon field at horizon crossing, we get an expression 


analogous to Eq. (59) of Ref. 


Pf°(k,l) 


20 (]: 


( 

1-84- \ ^ 

2 

Gfl\W 

w 


V 

z 0 1 

Tr (<5 + 


16 \ — <5 


3 w *?) r R (J^ + z)r R 


85 
97 2 


r R 


85 , 5 

9 ? + 2 


rS Pf\K 1 ) 


Ml) 


where 5 = ’ W = u+ ± l<T , Z = ^(l + Zcx — |), T r is the Gamma function and Gfg is 

a Meijer-G function. From this last expression we obtain the primordial power spectrum of the 
tachyon field. Some numerical values are shown in TABLE HI 
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